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I.  INTRODUCTION 


A.  BACKGROUND 

In  2015,  as  reported  in  the  Navy  Times  [1],  U.S.  Navy  Secretary  Ray  Mabus 
suggested  that  the  F-35  Joint  Strike  Fighter  “should  be  the  last  manned  strike  fighter 
aircraft  the  Department  of  Navy  will  ever  buy  or  fly.”  What  can  be  inferred  from  this 
statement  is  that  unmanned  aerial  vehicles  (UAVs)  will  continue  replacing  manned 
aircraft  in  essential  missions  and  in  support  of  ground  forces  to  accomplish  their 
missions.  Although  there  are  advantages  and  disadvantages  to  using  UAVs,  the  belief  that 
UAVs  will  eventually  replace  manned  vehicles  will  fuel  further  improvements  in  UAV 
technology.  One  of  these  areas  is  the  ability  to  charge  the  UAV  systems  without  having 
to  hardwire  a  connection  to  a  power  source.  To  achieve  this,  a  microwave  signal  from  a 
base  station  can  be  used  to  charge  the  UAV  once  it  lands,  as  shown  in  Figure  1.  Because 
of  the  close  proximity  of  the  UAV  to  the  base  station,  the  electromagnetic  beam  falls 
within  the  Fresnel  region.  This  wireless  power  transmission  (WPT)  approach  has  been 
studied  extensively  in  the  past  to  implement  this  concept  and  maximize  the  transfer 
efficiency. 


Figure  1.  Concept  Drawing  of  a  Base  Station  Array  Wirelessly  Transferring 

Power  to  a  Grounded  UAV 
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Wireless  power  transfer  has  been  explored  and  analyzed  using  both  an  inductive 
approach  and  a  radiative  approach  [2]-[9].  In  Figure  2,  the  block  diagram  for  the  WPT 
system  is  shown.  The  power  station  provides  the  power  that  is  used  to  charge  the  power 
plant  or  battery  in  the  client.  The  power  starts  off  as  low  frequency  A/C  power,  so 
conversion  to  radio  frequency  (RF)  is  needed  to  change  the  power  to  RF  so  that  the  signal 
can  be  transmitted  via  an  antenna  or  coil.  On  the  client  side,  the  rectifying  antenna 
receives  the  RF  power,  filters,  and  conditions  the  power  (transforms  its  voltage  and 
current  to  direct  current  (DC))  so  that  it  can  be  channeled  to  the  battery  or  power  plant. 


Base  Station 
(Master) 


Transmitting 


Recterma 


Client 

(Slave) 


Figure  2.  WPT  System  Block  Diagram  for  Battery  Charging. 

Source:  [2]. 


Wireless  power  transfer  has  gained  renewed  interest  since  the  advent  of  cell 
phone  technology  as  well  as  increased  drone  usage  in  the  military;  although  the  interest 
has  always  been  there,  the  minimization  of  circuit  components,  as  well  as  the  advent  of 
increased  processing  power  and  memory  size,  permits  the  implementation  of  efficient 
WPT  for  commercial  as  well  as  military  applications. 

From  a  smartphone  user’s  standpoint,  down  time  for  recharging  a  cell  phone 
means  the  user  has  less  time  to  conduct  business,  work  on  a  project,  or  just  leisurely  play 
a  game.  From  the  mission  standpoint,  if  the  operator  needs  to  recharge  the  drone,  then 
this  situation  equates  to  one  less  asset  to  fulfill  the  mission.  Any  downtime  can  have  dire 
consequences  for  ground  forces  that  rely  on  the  drones  to  be  their  eyes,  to  destroy  a 
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critical  enemy  outpost,  or  to  transport  equipment.  In  other  words,  the  longer  the  drone  is 
up  and  running,  the  more  we  can  support  troops  in  the  field. 

The  space  surrounding  an  antenna  can  be  broken  down  into  three  regions  [10]:  (1) 
reactive  near  field,  (2)  radiating  near  field,  and  (3)  far  field.  From  a  distance  of  0  to 

0.62^D3/A  ,  where  D  is  the  largest  dimension  of  the  antenna  and  A  is  the  wavelength, 

the  reactive  near  field  dominates.  From  0.62^  D3/ A  to  2  D2  j A  ,  “the  radiation  fields 
predominate  and  its  angular  field  distribution  is  dependent  upon  the  distance  from  the 
antenna.  The  field  pattern  is,  in  general,  a  function  of  the  radial  distance  and  the  radial 
field  component  may  be  appreciable”  [10].  This  region  may  not  exist  if  the  largest 
dimension  of  the  antenna  is  minimal  relative  to  the  wavelength.  This  region  is  also 
commonly  referred  to  as  the  Fresnel  region.  The  region  from  2D1  / A  to  infinity  is  known 
as  the  far-field  region.  In  this  region,  the  radiation  pattern  is  independent  of  distance  [11]. 
In  the  near  field,  because  of  the  dependency  of  the  radiation  pattern  on  the  distance  to  the 
antenna,  special  techniques  are  required  to  maximize  the  power  at  a  close  distance  to  the 
antenna. 

Technological  limitations  prevented  earlier  pioneers  in  the  WPT  realm  from 
implementing  efficient  concepts.  Since  rectification  is  a  nonlinear  process,  without 
precisely  knowing  the  power  level  at  a  given  distance  from  the  antenna,  it  is  very  difficult 
to  operate  the  rectenna  (a  special  type  of  antenna  that  is  used  to  convert  electromagnetic 
energy  into  DC  power  at  maximum  efficiency.  Lacking  an  efficient  method  for  delivery, 
WPT  is  not  practical.  Today,  technological  advances  enable  processors  to  compute  more 
quickly  and  begin  to  solve  these  efficiency  issues.  Commercial  software,  such  as  Savant, 
allows  the  user  to  simulate  the  radiated  field  distribution  in  the  near-field  environment  of 
the  antenna,  allowing  for  an  optimum  rectenna  design. 

B.  PREVIOUS  WORK 

Previous  work  has  been  done  in  the  analysis  of  near-field  focusing  for  WPT, 
specifically  in  [2].  The  WPT  method  was  introduced  in  the  1960s  for  vehicle  propulsion, 
and  researchers  from  Raytheon  eventually  developed  and  demonstrated  the  method  [12]. 
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With  the  successful  demonstration,  the  next  step  was  to  improve  the  efficiency  of  the 
power  delivered  as  well  as  the  form  factor  for  the  system.  Much  of  the  recent  research 
has  been  geared  towards  improving  these  two  criteria  in  the  near  field,  as  an 
improvement  in  processing  power  has  lent  itself  to  improvements  in  the  WPT 
performance. 

In  [3],  Ling  and  Yoon  surmised  a  new  analytical  upper  bound  for  power  transfer 
efficiency  (PTE)  for  WPT  using  a  pair  of  antennas  separated  by  a  spherical  material  shell 
and  showed  that  lossy  dielectric  materials  disturb  the  transverse  magnetic  (TM)  and 
transverse  electric  (TE)  mode  radiators.  Numerical  simulation  using  FEKO  validated  the 
theoretical  formulation,  where  FEKO  is  a  German  acronym  for  “FEldberechnung 
fiir  Korper  mit  beliebiger  Oberflache”  meaning  “field  calculations  involving  bodies  of 
arbitrary  shape.” 

Zhang  and  Cheng  [13]  experimented  with  two  efficient,  electrically  small  planar 
antennas  and  found  that  these  two  antennas  could  reach  high  PTE  without  adaptive 
matching  methods. 

In  the  commercial  world  of  cell  phones,  laptops,  tablets,  etc.,  Li  and 
Jandhyala  [14]  examined  the  use  of  a  retrodirective  antenna  array  for  WPT  with  the 
antenna  placed  at  a  static  location  charging  randomly  placed  devices  within  that  room. 
They  provided  design  guidelines  using  a  lookup  table  for  different  parameters:  element 
number,  element  size,  spacing,  and  focus  distance. 

Other  researchers  such  as  Deng  and  Kong  [4]  focused  their  efforts  on  improving 
the  components  in  the  WPT  system,  more  specifically  the  rectenna.  Using  a  common 
diode  they  simulated  the  rectenna  performance  using  their  new  method  and  theory  that 
provided  an  output  of  25.6  mW  and  52%  conversion  efficiency;  however,  researchers  in 
Japan  [6]  chose  to  develop  new  substrates  in  order  to  improve  the  RF-to-DC  conversion 
efficiency.  They  demonstrated  an  RF-to-DC  conversion  efficiency  of  74.4%  at  2.45  GHz 
at  a  5-W  input  power  using  a  10-finger  GaN  Schottky  diode.  Working  on  improving  the 
form  factor  of  the  rectennas,  researchers  in  Korea  [15]  designed  a  low-profile  printed 
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rectenna  in  the  X-band  that  showed  a  7%  and  4%  maximum  efficiency  for  a  two-by-two 
array  rectenna  and  a  four-by-four  rectenna,  respectively. 

The  researchers  from  Sendai  National  College  of  Technology  and  Tohoku 
University  [16]  calculated  the  transfer  efficiency  of  a  WPT  system  by  using  the  scattering 
parameters  (S-parameters)  of  the  WPT  system  modeled  as  a  two-port  lossy  network. 
Using  this  general  and  practical  method,  they  found  that  the  maximum  transfer  efficiency 
occurred  when  the  input  and  output  ports  are  both  conjugately  matched. 

On  the  same  note,  Koziel  and  Mongiardo  [17]  derived  a  surrogate-based 
optimization  method  for  improving  transfer  efficiency  using  equivalent  network  theory 
and  conjugate-image  impedances.  Although  a  full-wave  optimization  is  ideal  for  practical 
cases,  the  researchers  showed  that  the  surrogate  based  optimization  was  a  fair 
approximation. 

With  regard  to  the  near  field,  work  has  been  done  to  improve  the  transfer 
efficiency  using  several  adaptive  matching  methods,  specifically  simultaneous  conjugate 
matching,  frequency  tracking,  and  a  combination  of  the  two  [18].  Because  the 
simultaneous  conjugate  matching  method  was  not  practical  and  the  frequency  matching 
method  showed  poor  performance  beyond  the  strongly  coupled  region  of  the  two 
antennas,  the  hybrid  frequency  tracking  method  with  a  complex  load  matched  at  the 
target  distance  proved  to  have  stable  efficiency  across  the  near-field  region. 

At  the  Naval  Postgraduate  School,  research  in  the  area  of  WPT  has  gained 
momentum  over  the  years.  In  1998,  Suksong  [19]  discovered  that  a  monopole  on  a 
cylindrical  end  cap  provided  the  lowest  voltage  standing  wave  ratio  (VSWR)  and  a  nearly 
hemispherical  radiation  pattern  based  on  four  different  antenna  configurations.  He 
planned  to  use  this  antenna  to  remotely  power  a  micro  UAV  (MAV)  [19].  Continuing  on 
the  work  of  remotely  powering  the  MAV,  Vitale  [20]  studied  several  microwave  rectifier 
systems  and  eventually  achieved  a  maximum  RF-to-DC  conversion  efficiency  of  33%. 
He  also  showed  that  1.8-W,  1.3  GHz  microwave  signal  could  power  a  miniature  DC 
motor  in  free  space  at  a  distance  of  30  inches  from  the  transmitter  antenna.  In  2003, 
Tsolis  [21]  demonstrated  that  a  MAV  weighing  between  9  and  15  grams  and  deriving  the 
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components  from  previous  research  to  work  at  a  different  frequency  can  be  powered  via 
WPT  at  10  GHz.  Further  work  was  done  by  Tan  [22]  to  improve  on  several  subsystems 
(receiving  antenna,  pre-rectification  filtering,  rectifier,  and  post-rectification  filter)  of  the 
MAV  and  reduce  its  dimensions  and  weight.  He  realized  a  final  design  with  an  overall 
dimension  of  112.5  mnT,  an  estimated  weight  of  less  than  a  gram,  and  an  RF-to-DC 
conversion  efficiency  of  50%.  Toh  [23]  built  the  10  GHz  microstrip  rectenna  suggested 
by  previous  research  and  accomplished  a  conversion  efficiency  of  26%-37%  per  rectenna 
element,  a  fourfold  improvement  over  [21].  Adding  to  the  work  on  improving  rectenna 
design,  Liu  [24]  found  that  a  dipole  antenna  array  used  as  full-wave  rectenna  without 
low-pass  filter  showed  a  conversion  efficiency  of  65%  and  also  reduced  the  weight  of  the 
overall  system.  In  2012,  Huang  [25]  investigated  how  best  to  optimize  the  full-wave 
rectenna  and  obtained  a  maximum  efficiency  of  57%  at  10  GHz  with  a  200  mW  input. 

In  the  present  research,  we  build  on  the  previous  work  in  [2]  and  examine  several 
ideas  to  improve  the  efficiency  in  the  near  field  for  an  array-to-array  configuration. 

C.  OBJECTIVE  AND  APPROACH 

In  [2],  a  digital  phased  array  is  used  on  the  transmit  side  to  form  a  narrow  beam  to 
illuminate  the  client  receive  antenna.  There  are  many  caveats  in  order  to  optimize  the 
efficiency  in  transferring  power  to  the  client  antenna.  First,  when  a  receiving  point  is  in 
the  near-field  of  a  transmitter  antenna,  the  phase  for  each  element  in  the  array  needs  to  be 
corrected  to  create  a  focused  beam;  however,  the  correction  focusses  the  maximum  at  a 
point.  If  the  receiver  antenna  has  a  large  area,  or  is  very  close  to  the  transmitter  antenna, 
then  some  of  the  receiver  elements  will  be  outside  of  the  transmitter  beam.  Analysis  of 
the  client  antenna  must  determine  amplitude  and  phase  corrections  required  to  address 
this  issue  so  that  more  of  the  client  antenna’s  surface  area  can  be  coherently  illuminated 
by  the  focused  beam  to  ensure  more  power  can  be  transferred. 

One  proposed  solution  to  address  the  phase  error  problem  is  to  divide  the 
transmitter  array  into  subarrays  to  create  multiple  beams  that  can  be  individually  focused 
to  cover  more  of  the  client  antenna  aperture.  Creating  the  correct  excitations  to  synthesize 
each  subarray  beam  is  a  challenge  to  efficient  operation.  Amplitude  and  phase  corrections 


6 


on  the  client  antenna  for  each  subarray  beam  need  to  be  applied  so  that  all  beams  are 
combined  coherently.  All  these  requirements  can  be  simulated  using  Matlab  and  verified 
using  Savant  (a  commercial  antenna  modelling  software). 

In  this  thesis,  we  are  concerned  with  the  radiative  approach  to  the  charging 
problem  based  on  the  operating  range  (distance)  of  interest.  Typically,  when  the  client  is 
charging,  the  distance  from  the  charging  base  station  is  not  more  than  3.0  m,  so  the  far- 
field  condition  may  not  be  satisfied.  If  the  transmitter  array  uses  uniform  amplitude  and 
phase  weighting  across  its  elements,  the  full  gain  of  the  antennas  is  not  achieved  in  the 
near  field.  The  spherical  phase  error  correction  method  (also  referred  to  as  a  quadratic- 
phase  correction)  serves  to  solve  the  near-field  issue  of  unfocused  energy  in  the  Fresnel 
zone.  In  Figure  3,  a  typical  case  where  the  focused  array  achieves  a  12  dB  improvement 
over  the  unfocused  array  at  distances  less  than  3.0  m  for  a  transmitter  array  with  a  single 
receiver  dipole  at  300  MHz  is  shown.  In  this  thesis,  we  expand  the  radiative  approach 
beyond  a  single  receiver  dipole  by  considering  various  transmitter  and  receiver  array 
sizes.  Along  with  the  quadratic-focusing  method,  we  also  explore  the  focused  complex- 
conjugate  matched  method.  The  complex-conjugate  method  calculates  the  transmitted 
field  at  the  points  where  the  receiver  antenna  elements  are  located,  and  then  conjugate 
field  weights  are  applied  to  the  transmitter  antenna  elements. 
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Figure  3.  Savant  Simulations  for  Near-Field  Ranges  with  and  without 
Quadratic  Focusing.  The  Scenario  is  Power  Received  by  a  Single  Dipole 
with  a  Matched  Load  over  a  Ground  Plane  from  a  Ten-by-Ten  Array  with 
Ground  Plane  Transmitting  1.0  W  at  300  MHz.  Source:  [2]. 

The  approach  to  achieving  these  comparisons  is  first  to  develop  an  analytical 
model  for  the  antenna  transmission  so  that  parametric  dependencies  are  established.  After 
creating  the  model,  we  employ  simulations  using  computational  electromagnetic  (CEM) 
software  to  conduct  the  parametric  studies  and  determine  the  most  efficient,  or  optimum, 
configuration  and  operating  conditions. 

The  following  parameters  are  varied:  operating  frequency,  distance  between  the 
transmitter  and  receiver  arrays,  and  sizes  of  the  transmitter  and  receiver  arrays.  Several 
operating  frequencies  are  considered  to  examine  how  the  spherical  phase-error  correction 
method  affects  the  overall  power  received  by  the  UAV.  Using  the  same  dimensions  and 
extending  the  error  correction  to  include  conjugate  matching,  we  show  the  advantages  of 
using  conjugate  matching  over  a  simple  phase-error  correction.  Finally,  tradeoffs  are 
discussed  between  the  two  methods  to  determine  the  best  overall  method  for  WPT. 
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D.  THESIS  OUTLINE 


In  Chapter  II,  an  analytical  model  and  its  equations  are  developed  that  lay  the 
foundation  for  simulations.  The  spherical  error  correction  is  discussed,  as  well  as  the 
conjugate  matched  correction. 

In  Chapter  III,  we  present  the  approach  to  the  parametric  study,  a  Savant  software 
overview,  including  its  advantages  and  limitations,  and  the  simulation  process  and 
programming  required  to  automate  the  parametric  study.  The  simulation  results  for  the 
unfocussed,  quadratic  correction  and  conjugate  matched  cases  are  discussed. 

In  Chapter  IV,  the  results  are  summarized,  and  a  conclusion  for  the  study  is 
provided.  Future  areas  of  research  that  were  not  included  in  this  study  are  suggested. 
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II.  ANALYSIS 


A.  MODEL 

1.  System  Design 

Within  the  WPT  system,  the  master  station  transmits  a  continuous  wave  (CW)  to 
the  client.  For  the  master  station,  the  antennas  tend  to  be  of  the  high  gain  variety, 
focusing  continuously  at  the  client.  On  the  receiver  end,  the  UAV  antenna  typically  is 
low  gain  and  small  in  size  because  of  platform  limitations.  The  received  power  versus 
system  parameters  are  governed  by  the  Friis  equation  [26] 

p  _  P'G,GrA2L 
(4  Kdy 


where  Pt  is  the  transmitted  power,  Gr  is  the  charging  station  antenna  gain,  Gr  is  the 
receiver  (client)  antenna  gain,  and  A  =  c/  f  is  the  wavelength  ( c  is  the  speed  of  light,  / 
is  the  frequency  in  Hz).  The  factor  F  accounts  for  wave  propagation  effects  in  the 
medium,  such  as  multipath  and  attenuation,  L  is  an  efficiency  factor  that  includes  device 
losses,  rectifier  efficiency,  and  separation  d  is  the  center-to-center  spacing  between  the 
antennas. 


In  the  far  field,  the  antenna  gains  Gt  and  G,  are  given  by  the  general  formula  [26] 


G  = 


4^4 

/l2 


(2) 


where  Ae  is  effective  area. 


2.  Near  Field  Formulation 

a.  Transmitter  Array 

Typically,  when  the  UAV  is  charging,  the  receiver  antenna  on  the  UAV  and  the 
transmitter  antenna  on  the  charging  station  are  in  each  other’s  near  field.  If  no 
adjustments  are  made  to  the  antennas  and  uniform  amplitude  and  phase  weighting  is 

utilized,  then  the  full  gain  of  the  antennas  is  not  realized.  With  full  control  of  the  array 
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elements’  amplitude  and  phase,  the  charging  station  may  utilize  near-field  focusing.  Each 
element  can  be  fed  by  a  transmit  module  that  is  controlled  by  the  digital  beamforming 
processor.  The  processor  calculates  the  phase  difference  between  each  transmitter  array 
element  and  the  receiver  array  element  and,  thereby,  calculates  the  amplitude  and  phase 
weight  for  each  transmit  module.  By  doing  so,  the  transmitter  array  focuses  more  of  its 
energy  on  the  receiver  array  [2] . 

In  Figure  4,  the  near-field  geometry  is  shown.  For  a  general  case,  let  the  array 
focal  point  be  at  rf  [xf,  yf,zf )  and  the  observation  (field  point)  be  at  R(x,y,z)  .  The 

array  is  located  in  the  x-z  plane  and  has  a  rectangular  grid  with  Nx  by  N,  elements 

(1  <  m  <  Nx  and  1  <  n  <  N  _ )  spaced  dx  by  d, ,  respectively  [2].  The  normal  vector  to  the 
array  is  y . 

The  other  parameters  that  are  shown  in  Figure  4,  as  described  in  [2],  are  the 
position  vector  to  element  m,  n  ,  dmn  =  xmnx  +  zm„z ,  where 


mn  ^  x  9 

(3) 

.  _  2n-(N.  +  1)  , 

2  WZ’ 

and 

(4) 

y  =0. 

The  position  vector  to  the  focal  point  is  given  by 

(5) 

7f=xfx+yfy  +  zfz, 

and  the  position  vector  to  the  field  point  is 

(6) 

R  =  xx +yy  +  zz. 

The  vector  from  element  m,  n  to  the  focal  point  is 

(7) 
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(8) 


R,nn  =  (xf  -  xmn ) X  +  ( yf  -  ymn )  y  +  (zf  -  zmn ) I , 


and  the  vector  from  element  m,  n  to  the  field  point  is 

r„m  =  (x  -  xmn )  X  +  ( y  -  ymn )  y  +  (z  -  zmn )  z , 

where 


r 

mn 


(9) 


(10) 


Figure  4.  Planar  Array  Antenna  Near-Field  Model.  Source:  [2]. 


Let  the  elements  be  z  -directed  and  linearly  polarized.  A  general  form  for  the 
element  factor  is  [2] 
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Kn  (  0mn  ,  <f>,nn  )  =  ZC^  (  Tm  )  f  mn  (  #„,„  ,  </>mn  )  (11) 

where  fmn  {9mnAmn)  is  thc  normalized  element  pattern  as  a  function  of  angles  0mn , (f>mn , 
'i'(rmn)  =  e  lPr’""lrmn  >s  the  free  space  Green’s  function  with  P  =  2n  I X ,  C0  a  complex 
constant  that  depends  on  the  element  type  (e.g.,  dipole,  patch,  etc.),  0mn  is  the  angle 
between  the  element  m  ,  n  axis  (along  the  z  -axis)  and  the  field  point  position  vector 

cos#  =z-r  ,sin#  =  Jl  -  cos2  0  ,  (12) 

and  <f>mn  is  the  element  azimuth  angle  where 

^  (13) 

Xf~Xmn 

Normally,  all  elements  of  an  array  are  the  same,  exhibiting  the  relationships 
F  (9  ,<j>  )  =  F  (9  ,(h  )  and  f  (9  ,</>  )  =  f  (0  ,<j>  ).  If  the  elements  are  unit 

mn  \  mn’rmn/  \  mn’/mn/  J  mn  \  mn  ’  /  mn  /  J  \  mn’rmn/ 

current  amplitude,  z  -directed  half-wave  dipoles  at  height  It  above  a  ground  plane  that  is 
a  perfect  electrical  conductor  (PEC),  then  the  element  factor  is  [2] 

(n  Q  \ 

e-jf>rmn  cos  yC0S^J 

F  ( e„,n •  Kn )  =  ./1 201 - - -sin  (J3h sin  9nm  sin <f)mn ) .  (14) 

r  sin# 

mn  mn 

^  'J,,,,)  f(0mnA„n) 

The  ground  plane  factor  (the  last  sine  factor  in  Eq.  (14))  is  for  a  distance  in  the  far 
field  of  the  element.  This  factor  is  only  relevant  for  the  idealized  case  of  an  “isolated 
element  above  an  infinite  ground  plane.  The  ground  plane  factor  for  a  dipole  over  a  finite 
ground  plane,  in  an  array  environment,  or  in  the  near  field  will  differ  from 
sin  ( ph  sin  9mn  sin  (f>mn )  due  to  mutual  coupling  and  diffraction”  [2] . 

The  total  array  field  at  the  observation  point  is  given  by  the  weighted  sum 
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(15) 


E(R,04)  = 

m= 1  n= 1 

Nx  Nz  [  -iPr„m  ] 

=  ./120ZXX'1  - f  {emn^mn)Amnel<!>m"el¥m"  \ 

m= 1  n=l  [  r,im  J 

where  Amn  is  the  element’s  current  amplitude  coefficient,  is  the  phase  required  to 

focus  the  beam  at  r , ,  O . =  B  ( r,  -  r  )  ,  and  i// .  is  a  miscellaneous  calibration  and 

correction  phase. 

We  exploit  these  relationships  at  the  transmitter  as  well  as  the  receiver  to  improve 
the  received  power.  Next  we  describe  the  combination  of  signals  at  the  receiver. 

b.  Receiver  Array 

Let  a  receiver  array  be  centered  on  the  y  -axis  at  yp  with  numbers  of  elements  in 
the  x  and  z  directions  N'x  and  N[  and  spacings  d[  and  d[ ,  respectively  [2].  Similar  to 
the  transmitter  array,  the  element  locations  are  given  by  x'pq  =  [2p ~(N'X  +  l)d']/2  , 
z'pq  =  \_2q  -  (N[  + 1)  J']/ 2  and  y'pq  =0  (l  <  p  <  N'x ,  1  <  q  <  N[ )  as  shown  in  Figure  5. 

For  the  receiver  array,  we  use  primed  quantities  and  indices  to  differentiate  from 
the  transmitter  array  quantities  and  indices.  The  incident  field  at  element  p,q  is  obtained 
by  setting  its  location  as  the  field  point  in  Figure  4: 

E(R.e.^^E‘{x'lq,y'n,z'„)^En.  (16) 
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z  1 


Figure  5.  Total  Field  at  Element  p,  q  of  the  Receiver  Array.  Source:  [2]. 


Inserting  the  appropriate  values  into  Eq.  (15),  the  total  field  at  element  p,q 
affected  by  the  transmitter  array  elements  is  given  by 


Nx  Ny  f  ^-jPrmnpq 


E'Pq  =  i 1 20^E  E  I  - f  (  emnpq  >  <t>,nnPq  )  An 


(17) 


Y 

m= 1  n= 1  |  '  mnpq 


where 


rn,nPq  =  (XPq  -  Xmn  )  X  +  (  }' Pq  ~  J ,nn  )  9  +  (Zpg  ~  Zmn  ) 


Zpq  Zmn  I  z , 


(18) 


r  =\r  and  6  is  the  angle  of  the  axis  of  transmitter  element  m  ,  n  with  the  line 

mnpq  mnpqi  mnpq  &  ’ 

of  sight  (LOS)  vector  to  receive  element  p,q  (an  extension  of  Eq.  (12)),  and  [2] 

<19> 

^~mnPq 


By  the  same  token  and  using  Eq.  (13),  the  azimuth  angle  is 
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)=ZjlJA  <2°) 

X  —X 

pq  mn 

Receiver  element  p,  q  outputs  a  complex  voltage  given  by 

V  =  El  h  iff  ,</>'  )  (21) 

pq  pq  pq  y  mnpq  ’  T mnpq  J  v  ' 


where  hpq  is  the  effective  height  vector  of  receive  element  p,q  and  {d'mnpq,(/>'mnpq)  is  the 
direction  from  the  receiver  element  axis  with  the  LOS  to  transmitter  element  m  ,  n  with 


cos 


/  /V  \  _  ' mnpi 
y  wmnpq  J  I  _ 


V  *  7 

mnpq  pq 


and 


tan  (Cot)  =  -7 


r 

mnpq 


J  mn  y  pq 


X  —  X 
mn  pq 


(22) 


(23) 


The  relationship  between  the  receiver  and  transmitter  array  is  given  by 


A.#  A  A  »  A  A  »  A 

z  =Z,X  =-x,y  =-y 


(24) 


and  each  point  is  determined  from 

z  =  z,x'  =  -x,  y'  =  -y.  (25) 

In  the  Appendix,  the  derivation  of  the  effective  height  for  a  half-wave  dipole  on 
the  z!  -axis  above  a  PEC  ground  is  developed,  and  the  result  is 

<26> 

This  final  form  is  “an  idealized  case  of  an  isolated  element  above  an 
infinite  ground  plane”  [2].  For  a  “dipole  over  a  finite  ground  plane  or  in  an  array 
environment,”  [2]  mutual  coupling  and  diffraction  alter  the  effective  height.  Finally,  the 
“total  complex  voltage  at  the  output  of  the  receiving  array”  [2]  can  be  calculated  and 
determined  to  be  the  complex  sum  of  all  receiver  dipole  voltages 
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v  =  YYv  A'  ej¥'pq 

out  pq  pq 


N'  N'7  Nr  Ny 


=  j - YYYY\z-0'  h  iff  ,(/)'  )A'  ejy/'n 

J  q  I  pq  pq  \  mnpq  9  Y  mnpq  J  pq 

r*  p= 1  q= 1  m- 1  n= 1 

x'F  (r  )  /  ( 6  ,(b  )  A  ej<s>""eJ'/,"m  1 

y  mnpq  J  J  y  mnpq  9  Y mnpq  J  mn  I 


where  A!  and  w'  are  miscellaneous  calibration  and  correction  amplitudes  and  phases 

pq  '  pq  1  1 


applied  at  the  receiver,  and  0'pq  is  spherical  system  unit  vector  centered  at  p,q.  Given  that 
Eqs.  (23)  hold 

z-0'pq  =  -sm0'mnpq.  (28) 


‘For  both  arrays  with  half-wave  dipoles  above  a  PEC  ground  plane,”  [2]  the  final  result  is 


r*  p= 1  q= 1  I  m= 1  n= 1 


n 

\  COS  - W 

frmnpq  \  r\  mnpq 


\Phv,»nPq) 


n  , 

COS  ^WmnPq 

xA  ^  =^-sin  (Bhv'  ) 

mn  / 2  y '  mnpq  J 

Jl-(w') 

V  \  /  mnpq 


where  v  =sin#  sin  6  ,w  =cos  0  ,sin#  =,/l  -w2 

mnpq  mnpq  Y mnpq 9  mnpq  mnpp 9  mnpp  \l  mnpci  9 


mnpq 7  mnpq 


v'  =  sin6f  sin d>'  ,  and  w  =cos O'  ,sin#'  =  Jl  — (w') 

mnpq  mnpq  Y mnpq 9  mnpq  mnpq 9  mnpq  \l  \  /  mnpq 


As  can  be  seen  in  Eq.  (29),  a  summation  of  the  elements’  voltages  is  considered  for  series 
combination.  In  contrast,  when  the  dipoles  are  arranged  in  a  parallel  combination,  voltage 
remains  constant  and  results  in  Voul/(N[.N[)  .  The  time-average  power  is 
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(30) 


Keeping  only  the  first  two  terms,  we  get  a  quadratic  phase  error  of  A(.r) «  x1  jlR  . 

Substituting  the  quadratic  phase  error  into  Eq.  (15)  gives  e~lPx  '2R  as  the  phase  term 
because  the  R  +  A  path  is  longer  than  R\  therefore,  the  quadratic -phase  error  correction 

needed  is  elflx  2R .  Since  the  quadratic-phase  error  correction  only  requires  knowledge  of 
the  receiver  array’s  center,  this  method  requires  less  processing  than  conjugate-field 
matching. 
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When  the  distance  between  the  transmitter  and  receiver  arrays  is  large  (the  far- 
field  case),  the  spherical  wave  fronts  are  essentially  plane  waves,  and  the  phase  error  is 
minimal  across  the  array.  As  the  distance  is  shortened  between  the  two  arrays,  the  phase 
error  increases.  The  phase  error  correction  holds  as  long  as  the  transmitter  array  appears 
to  be  a  point  source;  however,  as  the  distance  between  the  transmitter  and  receiver  arrays 
shorten  even  more,  then  the  point  source  assumption  no  longer  holds  and  phase  error 
corrections  fail  to  improve  the  received  power.  The  complex-conjugate  method  provides 
the  actual  correction  needed  at  any  distance. 


Figure  6.  Path  Length  Difference  for  a  Point  Source  Illuminating  an 

Aperture 
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4. 


Conjugate-Matched  Correction 


In  [27],  it  was  shown  that  a  complex  conjugate  match  of  the  incident  focused  field 
resulted  in  an  increased  power  transfer  between  antennas.  Referring  to  Figure  5,  we 
consider  a  receiver  element  ( p,q )  at  distance  y  .  The  transmitter  array  has  excitation 

coefficients  Gmn  ,  and  the  receiver  array  has  coefficients  A  .  The  transmission 
coefficient  between  element  (, m,n )  and  (p,q)  is 


8„ 


_  V Asinfl, 


mnpq^r  ^ \nnpq 


^mnpq^~ 


(35) 


where  A(  and  Ar  are  the  effective  areas  of  each  arrays’  elements.  We  desire  the  Gmn  for 
the  field  that  is  conjugate  matched  at  the  receiver  plane  ypo  : 

G*  =VVg  A*  e~iPrmnpq.  (36) 

mn  o  mnpq  pq  v  / 

P  q 

We  apply  these  weights  to  the  transmitter  antenna.  Extra  processing  is  required,  and  it  is 
necessary  to  know  the  location  of  the  receiver  array  elements  relative  to  the  transmitter 
array  elements.  The  complex-conjugate  match  method  provides  the  correction  for  a 
specified  geometry  (spacing  between  arrays,  array  parameters,  frequency,  etc.).  As 
opposed  to  quadratic-phase  correction,  the  conjugate  method  utilizes  both  amplitude  and 
phase  corrections. 
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III.  SIMULATION 


A.  INTRODUCTION 

1.  Parametric  Study  Approach 

The  main  goal  of  this  research  is  to  increase  the  power  received  by  the  receiver 
array  by  correcting  the  phase  of  the  transmitter  array  or  by  correcting  both  the  phase  and 
amplitude  of  the  transmitter  array.  Part  of  the  parametric  study  revolves  around 
confirming  the  results  from  [2]  by  varying  the  distance  between  a  ten-by-ten  transmitter 
array  and  a  single  receiver  dipole.  Once  this  is  accomplished,  an  array  replaces  the 
dipole  at  the  receiver  end. 

With  arrays  at  both  the  transmitter  and  receiver,  we  vary  parameters  to  determine 
which  configuration  increases  the  power  received.  Since  the  process  was  automated 
(discussed  below),  the  uniform,  spherical-phase  correction  method,  and  conjugate-match 
method  are  plotted  on  the  same  graph.  Distance  is  varied  from  zero  meters  to  two  meters 
since  anything  beyond  two  meters  is  considered  in  the  far  field  for  our  frequencies  of 
consideration. 

The  size  of  each  of  the  arrays  are  parameters  of  consideration  for  the  effect  on  the 
power  reception.  One  of  the  array  sizes  is  kept  constant  while  the  other  array  size  is 
changed.  For  example,  if  the  size  of  the  transmitter  array  is  varied,  then  the  receiver  array 
has  a  constant  size.  The  transmitter  array  was  varied  to  be  both  smaller  and  larger  than 
the  receiver  array.  In  our  research  the  arrays  both  were  kept  square  for  simplicity.  Due  to 
the  limitation  of  the  available  PC,  we  only  investigated  arrays  up  to  30-by-30. 

Finally,  the  frequency  was  adjusted  for  each  of  the  scenarios  discussed  above.  We 
focused  our  research  primarily  in  the  Industrial,  Scientific,  and  Medical  (ISM)  bands  and 
the  X-band  (specifically  2.4  GHz,  5.8  GHz,  and  10  GHz).  The  ISM  band  encompasses 
most  of  the  commercial  applications,  and  the  X-band  range  covers  radar  applications. 
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2. 


Savant 


Savant  is  an  antenna  simulation  software  that  takes  into  account  the  installation 
environment  such  as  a  tank,  car,  building,  or  other  platform.  It  can  predict  the  radiated 
field  distributions  in  the  environment  surrounding  an  antenna,  the  far-field  radiation  and 
receive  patterns  of  antennas,  and  the  cosite  coupling/isolation  between  antennas  mounted 
on  the  same  platform  or  two  separate  platforms. 

Savant’s  simulation  method  is  based  on  the  shooting-and-bouncing  ray  (SBR) 
method.  Savant’s  SBR  is  a  combination  of  two  methods:  geometrical  optics  (GO)  and 
physical  optics  (PO).  GO  describes  the  simple  ray-tracing  techniques  that  have 
traditionally  been  used  in  optics.  In  order  to  use  GO,  it  is  assumed  that  wavefronts  are 
locally  planar  and  waves  are  TEM  [28]. 

PO  describes  the  approximation  of  the  current  induced  “using  ray  optics  to 
estimate  the  field  on  a  surface  and  integrating  that  field  over  said  surface  to  calculate  the 
transmitted  and  scattered  field”  [29] .  The  PO  induced  current  is  [30] 

y  =J2 nsxH‘,  lit  region 
*  |  0,  shadow 

where  ns  is  the  normal  to  the  surface.  As  shown  in  Figure  7,  H'  is  the  magnetic  field 

incident  on  the  scatterer.  Using  the  radiation  integral  with  the  PO  induced  current,  we 
obtain  the  scattered  field, 

Es  (r)  =  -22120  [ [  k  x  k  x  \2ns  x  H'  1  dS  .  (38) 

An  J  r  —  r'\ 

s  I  ! 

The  position  vectors  to  the  source  point  r  and  observation  point  r  are  shown  in 
Figure  7. 
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Observer 


Savant  launches  the  GO  rays  from  the  transmitter  antenna  towards  the  surface  and 
determines  which  areas  directly  illuminated.  The  areas  that  are  illuminated  induce  the  PO 
currents  that  “radiate  to  far-field  observation  angles,  near-field  observation  points,  or 
receive  antennas  to  generate  scattered  field  contribution”  [31].  These  reflected  GO  rays 
travel  and  hit  other  surfaces  in  the  installed  environment,  inducing  secondary  PO 
currents.  As  this  process  continues  onto  the  next  point  and  the  next,  Savant  calculates 
multiple  bounce  scattering  and  determines  the  scattered  fields  for  arbitrary  shapes. 

For  antenna  coupling,  the  expression 


b  =  ,jA  E  ■  d„ 


(39) 


is  used  to  determine  the  receive  antenna’s  coupled  power,  where  dRx is  the  scaled  far- 
field  directivity  vector  of  the  receive  antenna,  E  is  the  incident  or  scattered  electric  field 
at  the  receiver  antenna,  and  r/Q  is  the  free-space  impedance.  The  quantity  b  is  the 
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complex  amplitude  of  the  received  signal  delivered  to  a  transmission  line  that  is 
impedance  matched  to  the  receive  antenna. 

Savant  provides  the  user  with  easy  to  use  functions.  It  accepts  multiple  computer- 
aided  design  (CAD)  file  formats,  so  whatever  format  used  in  the  platform  design  can  be 
accommodated  and  imported  into  Savant  for  simulation.  Savant  also  comes  with 
preloaded  antenna  shapes  and  sizes  so  a  new  user  can  run  simulations  on  many  common 
antennas.  If  the  user  has  an  antenna  described  in  a  Savant  antenna  element  list  (SARR) 
file,  then  Savant  uses  that  file  to  generate  its  simulations.  The  three  analysis  types  provide 
an  easy  to  use  interface  for  beginning  the  analysis.  The  user  can  create  far-field  patterns 
and  field  distributions  (  E  and  H  fields)  for  installed  transmit  antennas.  Coupling 
coefficients  of  the  installed  transmitter  and  receiver  antennas  can  be  derived.  Once  the 
analysis  is  completed,  the  color-coded  visualization  of  the  fields  provides  a  tool  for 
understanding  where  the  power  is  distributed  around  the  platform  and  environment. 

Savant  allows  their  application  to  be  integrated  into  Matlab.  The  developers 
provided  an  application  programming  interface  (API)  for  Matlab  to  run  Savant  via 
Matlab  code.  This  functionality  allowed  for  the  parametric  study  to  be  automated.  Savant 
also  provided  feedback  on  each  run  made,  listing  the  details  of  its  simulation. 

Although  Savant  is  relatively  fast  for  analyzing  antennas,  the  computational 
capability  of  a  typical  Windows  desktop  had  a  limit  of  20-by-20  element  arrays  for  the 
transmitter  and  receiver  antennas.  For  a  computer  with  a  3.4  GHz  i7  Intel  quad  core 
processor  and  24  GB  of  memory,  each  configuration  can  take  up  to  six  hours  to  complete, 
limiting  the  scope  of  the  research.  Because  Savant  uses  SBR  as  its  basis  for  analyzing 
antennas  in  their  surroundings,  the  results  are  still  an  approximation  of  the  real  world. 
Other  simulation  techniques  and  measurements  of  the  transmitter  arrays  and  receiver 
arrays  may  not  exactly  match  the  results  of  these  simulations. 

3.  Near  Field  vs.  Far  Field 

In  [2],  Jenn  presented  the  case  that  the  focused  array  transmitter  antenna  to  a 
single  dipole  receive  antenna  showed  an  improvement  over  the  unfocused  array 
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transmitter  antenna.  As  part  of  this  thesis  we  reaffirmed  the  results  from  [2]  and  present  it 
in  Figure  8.  Although  it  is  a  simple  case,  the  results  shown  in  Figure  8  demonstrate  the 
improvement  possible  using  corrections.  An  operating  frequency  of  1  GHz  was  used  for 
the  simulation.  The  transmitter-array  element  spacing  is  half  a  wavelength  apart  and 
quarter  of  a  wavelength  above  the  ground  plane.  The  receiver  dipole  antenna  is  half  a 
wavelength  long  and  is  also  a  quarter  of  a  wavelength  above  its  ground  plane.  The 
scenario  is  depicted  in  Figure  9,  with  the  transmitter  array  and  dipole  1.0  m  apart.  As  can 
be  seen  in  Figure  8,  the  phase  correction  shows  an  improvement  for  the  Fresnel  field 
range.  Within  the  1.0  m  range,  the  focused  transmitter  array  delivers  almost  15  dB  more 
power  to  the  receiver  dipole  as  opposed  to  the  unfocused  array.  Beyond  3.0  m  there  is  no 
significant  difference  between  the  focused  and  unfocused  array  because  the  far-field 
conditions  begin  to  take  hold.  Although  the  phase  differences  cause  the  focused  and 
unfocused  arrays  to  have  different  patterns,  the  radiation  pattern  shape  in  the  far  field  is 
independent  of  distance.  This  case  provides  a  good  indicator  for  the  improvement  that 
can  be  achieved  using  a  focused  transmitter  array. 

From  Figure  10,  we  present  the  more  general  Savant  scenario  of  two  arrays.  The 
transmitter  and  receiver  array  elements  are  spaced  is  half  a  wavelength,  and  the  elements 
are  above  the  ground  plane  by  a  quarter  of  a  wavelength.  With  Matlab  automation 
discussed  below,  the  operating  frequency  and  distance  between  the  transmitter  and 
receiver  array  varies  depending  on  the  Matlab  created  Savant  configuration  file.  The 
ground  planes  for  both  arrays  extend  a  quarter  of  a  wavelength  beyond  the  array 
elements.  The  transmitter  array  elements  have  weights  that  are  dependent  on  the  SARR 
file.  On  the  other  hand,  the  receiver  array  elements  use  the  default  (uniform)  weights 
from  Savant. 
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Figure  8.  Coupling  of  a  Ten-by-Ten  Transmitter  Array  to  a  Single 
Dipole  Receiver  Antenna  at  1  GHz 


Figure  9.  Savant  Scenario  of  Ten-by-Ten  Transmitter  Array  to  Single 

Dipole  at  1  GHz 
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The  transmitter  array  takes  advantage  of  Eqs.  (3)-(5)  and  adjusts  the  phase  from 
each  element  to  the  observation  point  along  the  y-axis  in  Figure  4.  By  doing  this  the 
transmitter  array’s  main  beam  is  focused  at  the  dipole  (as  opposed  to  infinity  for  uniform 
phase)  and  a  majority  of  the  energy  is  incident  on  the  receiver  dipole. 


Figure  10.  Savant  Scenario  of  Transmitter  Array  Weighted  SARR  File 
and  Receiver  Array  Using  Default  Savant  Weights 

4.  Matlab  Automation 

In  simulating  a  transmitter  array  to  a  single-dipole  scenario,  there  was  code  to 
work  with  from  [2].  The  code  creates  the  SARR  files,  which  are  Savant  array  element  list 
files  with  weights  of  the  transmitter  array  in  the  Savant  scenario.  The  weights  are 
determined  by  the  frequency,  number  of  elements,  and  distance  between  the  transmitter 
array  and  the  receiver  array.  With  these  parameters,  the  phase  corrections  are  calculated 
and  stored  in  the  SARR  file.  For  focusing  using  quadratic  corrections,  only  phases  are 
required.  For  conjugate  corrections  both  amplitude  and  phase  are  required. 
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Savant  allows  the  user  to  use  preloaded  antennas  or  antennas  from  file,  so  the 
transmitter  array  utilizes  the  “antenna  from  file”  feature  by  loading  the  SARR  file  to 
represent  the  elements  of  the  array.  The  receiver  array  consists  of  the  preloaded 
rectangular  array  with  a  ground  plane.  Once  loaded  Savant  creates  the  transmitter  array 
from  the  SARR  file  and  creates  the  receiver  array  from  its  library  as  shown  in  Figure  10 
and  Figure  11.  Although  Figure  9  and  Figure  10  appear  similar,  the  transmitter  antenna  is 
an  array  and  a  dipole,  respectively.  Also,  the  green  plates  are  the  ground  planes  for  the 
two  arrays  and  are  set  to  be  quarter  of  a  wavelength  behind  the  arrays.  In  Figure  11,  the 
distance  between  the  transmitter  and  receiver  array  configurations  is  0.5  m  and  is  shown 
for  clarity.  Savant  also  requires  the  type  of  antenna  representation  for  each  element  of  the 
array,  so  for  this  case,  current  sources  were  chosen  for  the  dipole  elements,  and  we  only 
considered  one  reflection  off  the  ground  plane. 


Figure  11.  Savant  Scenario  of  a  Three-by-Three  Transmitter  Array  to 
a  Three-by-Three  Receiver  Array 

The  parameters  were  changed  automatically  to  run  each  scenario.  Parameters  that 
must  be  changed  include  the  distance,  number  of  elements  and  frequency.  For  each  case  a 
new  SARR  file  was  created  to  focus  the  beam  correctly;  moreover,  each  change  in  array 
size  demanded  a  change  in  the  ground  plate  dimensions.  This  type  of  effort  was  required 
for  both  the  unfocused  transmitter  array  as  well  as  the  focused  transmitter  array. 
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Depending  on  the  size  of  the  arrays,  each  run  required  five  to  30  minutes  of  computing 
time  using  the  computer  stated  earlier  for  each  datapoint.  Automation  saves  time  in 
running  a  large  number  of  scenarios. 

Savant’s  scenarios  are  stored  in  a  configuration  file  based  on  extensible  markup 
language  (XML).  XML  is  a  markup  language  like  the  hypertext  markup  language 
(HTML)  with  an  emphasis  on  storing  and  transporting  data  as  opposed  to  displaying  data. 
In  the  case  of  Savant,  the  configuration  file  contains  tags  for  defining  various  aspects  of 
the  scenario  such  as  antenna  type,  length  of  the  array,  frequency,  etc.  These  different 
configuration  tags  are  organized  into  nodes  and  tagged  as  “<NODE_CREATE>.”  For 
each  antenna  and  ground  plate  node,  there  are  dimension  tags  that  exist  to  describe  the 
length,  width,  and  position  of  the  array.  In  addition,  the  spacing  between  the  elements  of 
the  array  is  modifiable.  In  order  to  change  the  SARR  file,  the  text  within  the  transmit 
array  tag  “<Filename>“  was  modified  for  the  unfocused  or  focused  SARR  file.  Lastly,  to 
change  the  frequency,  the  data  stored  in  the  following  tags  were  modified:  <Frequency>, 
<ResonantFrequency>,  “DesignFrequency>,  <SingleFrequency>,  and  <MaxFrequency>. 
All  modifications  were  done  via  Matlab  since  Matlab  provides  functions  to  create, 
modify,  and  delete  XML  files. 

There  are  many  parameters  that  can  be  adjusted  to  examine  which  method 
provides  the  best  power  reception.  In  this  thesis  the  focus  was  on  the  frequency,  distance 
between  transmitter  and  receiver  antennas,  and  the  number  of  array  elements.  The 
spacing  between  elements  in  both  arrays  was  kept  to  half  a  wavelength,  so  the  frequency 
determined  the  spacing.  This  is  an  area  that  can  be  explored,  but  due  to  time  limitations, 
we  chose  not  to  focus  attention  on  spacing.  Because  of  the  limitations  of  processing 
power  of  the  available  computer,  array  sizes  were  kept  below  30-by-30.  As  mentioned 
before,  even  with  a  30-by-30  array,  the  Savant  simulation  took  six  hours  to  complete 
because  of  the  increasing  time  to  calculate  bounces  for  every  additional  element. 
Additionally,  the  shapes  of  the  arrays  were  restricted  to  a  square  configuration.  We 
choose  three  different  operating  frequencies:  2.4  GHz,  5.8  GHz,  and  10  GHz.  The  first 
two  frequencies  are  based  on  the  ISM  band.  ISM  band  consists  of  many  common 
household  items  such  as  cellular  signals,  Bluetooth,  and  Wi-Fi.  Air-to-air  applications 
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such  as  radar  warning  receivers  (RWR)  and  jammers  commonly  operate  in  the  latter 
frequency  band.  Simulations  were  done  by  varying  the  distance  between  the  arrays  from 
0.1  m  to  2.0  m.  These  distances  illustrated  the  effectiveness  of  the  methods  from  near¬ 
field  to  far-field.  Beyond  2.0  m,  the  far-field  behavior  exhibits  inverse  distance- squared 
behavior. 

B.  SIMULATION  RESULTS 

The  results  shown  in  this  section  are  for  several  scenarios  run  on  Savant  using 
unfocused,  focused,  and  complex-conjugate  element  excitations  on  the  transmitter 
antenna. 

1.  Transmitter  Array  to  Receiver  Array 

Since  the  receiver  array  is  limited  by  the  size  of  the  UAV,  the  maximum  number 
of  elements  is  influenced  by  the  choice  of  the  frequency  of  operation.  For  small  UAVs, 
typically  the  room  allowed  for  the  receiver  array  is  about  6.0  inches;  therefore,  for 
10  GHz,  the  largest  receiver  array  size  is  about  ten-by-ten  elements.  For  the  5.8  GHz 
operating  frequency,  the  largest  array  is  six-by-six,  and  for  the  2.4  GHz  frequency,  the 
largest  array  is  two-by-two. 

Two  different  methods  are  presented  and  compared  to  determine  their  effect  on 
the  power  received.  As  was  done  with  the  single  dipole  receiver  antenna,  the  quadratic- 
phase  correction  approach  is  the  first  method  examined  for  its  effectiveness.  A  further 
improvement  on  the  focused  approach  is  to  take  the  complex  conjugate  of  the  field  at  the 
receiver  element  locations,  thereby  providing  a  complete  match  and  delivering  more 
power  to  the  receiver  array.  The  complex-conjugate  method  was  proven  very  effective 
in  [27].  The  uniform  method  is  used  as  a  baseline  to  compare  the  effectiveness  of  the 
correction  methods.  As  shown  below,  Figures  12  through  20  are  plots  of  the  antenna-to- 
antenna  coupling  Pr/Pt  when  Pt  =1.0  W.  The  case  shown  in  Figure  12  verifies  that  all 
three  excitations  give  the  same  result  in  the  far  field.  For  two-by-two  arrays  at  2.45  GHz, 
the  far-field  condition  is  satisfied  at  about  2D2/T  =  2(0.0625)‘/0.125  =  0.0625  m. 
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Figure  12.  Power  Received  for  a  Transmitter  Array  (Two-by-Two) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Two-by-Two)  at  2.4  GHz 


2.  Simulation  of  the  Correction  Methods 

The  focused  spherical  (quadratic)  correction  method  requires  only  knowledge  of 
the  location  of  the  center  of  the  receiver  array.  This  simplifies  the  processing  and 
hardware  at  the  transmitter  array.  Based  on  the  study,  the  size  of  the  transmitter  array 
needs  to  be  at  least  twice  the  size  of  the  receiver  array  for  the  quadratic-correction 
method  to  be  effective.  As  can  be  seen  from  Figures  13  and  14,  the  quadratic-phase 
correction  method  does  not  work  well  if  the  transmitter  array  is  the  same  size  or  is 
smaller  than  the  receiver  array.  From  Figures  15  through  17,  when  the  transmitter  array 
dimensions  are  at  least  a  factor  of  three  greater  than  the  receiver  array,  the  power 
delivered  at  0.5  m  and  more  is  increased  by  3  dB  or  greater. 

The  quadratic-correction  method  does  not  work  well  when  the  distance  between 

the  transmitter  array  and  receiver  array  is  too  small.  The  steep  drop  off  in  power  received 

can  be  seen  in  Figures  13  through  17.  From  the  receiver  array  observation  point,  the 
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transmitter  array  no  longer  looks  like  a  point  source,  so  the  quadratic-phase  correction 
method  fails  to  improve  the  received  power.  The  junction  point  can  be  seen  where  the 
power  received  in  the  uniform  and  quadratic-phase  correction  methods  are  equal  in 
Figures  16  and  17.  Once  past  this  junction  point,  however,  the  quadratic -phase  correction 
and  conjugate-match  method  both  improve  the  received  power.  As  can  be  seen  from 
Figure  16,  from  1.4  m  and  on,  both  methods  are  equally  better  than  uniform  excitation. 

The  complex-conjugate  method  has  the  best  results  in  regards  to  power  delivered 
to  the  receiver  array.  The  power  delivered  does  not  drop  off  when  the  distance  between 
the  transmitter  array  and  the  receiver  array  is  too  small.  Even  when  the  transmitter  array’s 
dimensions  are  equal  or  less  than  that  of  the  receiver  array,  the  conjugate-excited 
transmitter  array  maintains  at  least  0.5  dB  advantage  over  the  uniform  method. 


Distance.(m) 

Figure  13.  Power  Received  for  a  Transmitter  Array  (Six-by-Six) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Six-by-Six)  at  5.8  GHz 
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Figure  14.  Power  Received  for  a  Transmitter  Array  (Ten-by-Ten) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a 
Receiver  Array  (Ten-by-Ten)  at  10  GHz 


Distance, (m) 

Figure  15.  Power  Received  for  a  Transmitter  Array  (Six-by-Six) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Two-by-Two)  at  2.4  GHz 
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Figure  16.  Power  Received  for  a  Transmitter  Array  (24-by-24) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Six-by-Six)  at  5.8  GHz 
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Figure  17.  Power  Received  for  a  Transmitter  Array  (30-by-30) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Five-by-Five)  at  10  GHz 
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Figure  18.  Power  Received  for  a  Transmitter  Array  (15-by-15) 

Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Two-by-Two)  at  2.4  GHz 

The  conjugate  method  requires  knowledge  of  where  all  the  elements  are  in  both 
transmitter  and  receiver  arrays.  The  extra  processing  power  and  hardware  needed  may 
not  be  worth  the  extra  power  delivered  for  distances  less  than  0.5  m;  however,  if  the 
receiver  array  is  stationary,  then  this  only  needs  to  be  determined  once  and  is  up  to  the 
sponsors  of  the  UAV  program  to  decide.  This  is  even  more  evident  when  we  look  at 
Figures  15  through  17.  There  is  at  least  an  8.0  dB  advantage  over  the  uniform  method  for 
small  distances. 

The  improvement  in  power  over  distance  is  proportional  to  the  ratio  of  the 
transmitter  array  to  the  receiver  array  Nt/Nr ,  where  Nt  is  the  length  of  the  transmitter 

array  and  Nr  the  length  of  the  receiver  array.  For  example,  if  we  look  at  Figure  15  and 

Figure  18  at  1.0  m,  the  power  received  is  -6  dB  and  -1  dB,  respectively.  If  the  UAV  is 
1.0  m  away  from  the  charging  base  station  operating  at  2.4  GHz,  then  we  only  need  to 
design  a  20-by-20  transmitter  array  using  quadratic-phase  correction.  If  the  UAV  requires 
charging  at  a  much  closer  distance,  such  as  0.2  m,  then  we  would  design  a  four-by-four 
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transmitter  array  using  quadratic-phase  correction.  The  conjugate-match  method  provides 
the  most  versatile  solution  since  the  quadratic-phase  correction  method  fails  at  close 
range  when  the  transmitter  array  no  longer  looks  like  a  point  source.  This  research  gives 
the  designer  guidelines  to  determine  what  method  of  correction  is  advantageous  given 
their  transmitter  and  receiver  array  sizes. 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2 

Distance, <jn) 

Figure  19.  Power  Received  for  a  Transmitter  Array  (20-by-20) 
Excited  by  Uniform,  Quadratic,  and  Conjugate  Coefficients  at  a  Receiver 
Array  (Two-by-Two)  at  2.4  GHz 


3.  Comparison  of  Analytical  Formulas  and  Simulation 

A  simulation  was  performed  to  compare  the  analytical  results  obtained  using  Eq. 
(31)  to  simulated  results.  The  data  in  Figure  20  compares  the  power  received  using  the 
focused  (quadratic  correction)  and  unfocused  cases.  As  can  be  seen  in  the  near  field 
range,  the  focused  case  provides  an  improvement  over  the  uniform  (unfocused)  case.  The 
curves  fall  in  line  with  the  behavior  observed  in  the  other  simulations.  As  long  as  the 
transmitter  array  is  seen  as  a  point  source,  either  focused  methods  may  be  used  to 
improve  power  reception. 
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Figure  20.  Power  Received  by  a  Single  Dipole  over  a  Ground  Plane 
from  a  Ten-by-Ten  Array  with  Ground  Plane  Transmitter  1  W  at  300  MHz. 
Comparison  of  Summation  Formula  to  Savant  Simulations  for  Near-Field 
Ranges  with  and  without  Focusing.  Source:  [2]. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Our  research  is  an  extension  of  the  technique  proposed  in  [2]  by  considering  two 
methods  for  correcting  the  phase  errors  between  elements  of  the  transmitter  and  receiver 
arrays  in  the  near  field.  Using  Savant  and  Matlab,  we  conducted  a  parametric  study  to 
determine  which  of  the  two  methods  increased  the  power  received  the  most.  By  varying 
the  sizes  of  the  transmitter  and  receiver  arrays,  changing  the  frequency  of  operation,  and 
varying  the  distance  between  the  arrays,  we  discovered  that  the  conjugate-match  method 
consistently  outperformed  the  spherical-phase  correction  method.  Using  the  conjugate- 
match  method,  we  found  that  the  power  received  by  the  receiver  array  increased  by  as 
much  as  12.0  Db,  as  long  as  there  are  a  sufficient  number  of  transmitter  array  elements  as 
compared  to  receiver  elements.  Although  the  conjugate  method  is  far  superior  to  the 
spherical  phase-correction  method,  it  requires  knowledge  of  the  all  the  receiver  element 
locations  as  well  as  calculation  of  the  amplitude  and  phase  weights.  All  these 
requirements  mean  more  complex  hardware  and  processing  power. 

If  the  receiver  array  is  stationary,  then  the  processor  only  needs  to  calculate  the 
locations  of  the  receiver  elements  once.  Since  the  UAV  is  more  likely  than  not  stationary 
during  charging,  the  conjugate-match  method  is  a  viable  method  for  improved  power 
transmission.  With  the  ever  increasing  improvement  in  processing  power,  calculating  a 
moving  UAV  array’s  receiver  element  locations  may  prove  trivial  enough  that  the 
conjugate-matching  method  can  be  used  transfer  power  to  a  moving  client. 

Due  to  time  constraints,  several  parameters  were  not  investigated.  In  our 
scenarios,  the  transmitter  and  receiver  arrays  lined  up  center  to  center  with  no  offset.  This 
is  an  area  that  should  be  researched.  Splitting  the  transmitter  array  into  sub-arrays  and 
focusing  each  sub-array  on  the  center  of  the  receiver  array  is  the  next  step.  As  mentioned 
earlier,  the  spacing  between  the  elements  for  the  arrays  was  kept  to  a  half  a  wavelength 
and  exploring  the  effect  of  other  spacings  is  worth  pursuing.  Changing  the  antenna  type 
from  a  dipole  is  another  study  that  could  prove  useful.  Also,  varying  the  length  and  width 
of  the  arrays  rather  than  having  square  arrays  might  improve  the  power  reception.  A 

change  in  the  CEM  software  is  worth  investigating.  Full-wave  solvers  such  as  the  method 
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of  moments  can  include  the  mutual  coupling  between  the  array  elements.  The  GO  ray 
tracing  does  not  rigorously  model  the  mutual  coupling  between  the  arrays,  which  can  be 
important  at  small  distances.  Also,  cloud  computing  should  be  utilized  to  speed  up  the 
simulations. 
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APPENDIX 


Consider  element  m,n  to  be  a  z  -directed  dipole  at  a  height  h  above  a  ground 
plane  that  is  located  in  the  x-  z  plane  [2] .  The  ground  plane  factor  is  [2] 


GPF  =  eJ 


2./sin(/?/r  cos  ay  j 


(40) 


where  cos  a  =v  =  sin  0  6  is  the  direction  cosine.  The  element  factor  for  a  z  -dipole 

ymn  mn  mnrmn  r 

over  the  ground  plane  is  [2] 


EF  =  z2  j  sin  ( Bh  sin  0nm  sin  d>  ) 

mn  •’  \  >  mn  T  mn  / 


-}Prm 
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=  -120z- 


sin  {phvmn)- 
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-  vv„ 


The  effective  height  is  determined  from 


V„„,  =  E  h 


(42) 


and 


E  =  -  jn—^—e~JI}rh  ,  (43) 

mn  J  I  a  mn  ’  '  ' 

47rr 

where  hmn  is  the  effective  height.  For  a  half-wave  dipole  with  a  ground  plane,  the 
effective  height  is 
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\/  mn  t  mn )  mn 


(44) 


Typically,  we  choose/?  -XI 4  so  that  [3h  =  n  /  2 .  This  is  an  idealized  far-field  result  for 
an  isolated  element  above  an  infinite  ground  plane.  For  a  dipole  in  an  array  setting  above 
a  finite  ground  plane  or  in  the  near  field,  mutual  coupling  and  diffraction  affects  the 
effective  height  and  may  not  result  in  a  sin  (/?/?  sin  6mn  sin  </>nm )  factor.  An  efficiency 

factor  A  and  taper  exponent  B  can  be  added  to  the  ground  plane  factor  to  characterize  the 
behavior  in  the  environment: 

A sinB  sin  ( [ih  sin  6mn  sin  <f)mn ) .  (45) 
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